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Abstract — We derive a system of propagation 
equations in a Double Negative (DN) material in a way 
that appears to differ from previous derivations— 
although the end result is the same. Our derivation 
assumes the Poynting vector theorem applies, real 
materials always have some loss, 4œ) and /(@) are 
obtained from real materials, and wave energy traveling 
in a specified direction must always be accompanied by a 
loss of energy in that direction. Additional mathematics 
beyond Maxwell's equation is not required. Energy 
losses per unit length of travel are finite, and can be 
extremely small. Propagation in a lossless DN media is 
found as the mathematical limiting solution of an 
extremely small energy loss per unit length. When 


developed along these principles, the equations 
developed for designing leaky antennas are 
straightforward. 


Index Terms— Double Negative material, Poynting 
vector theorem, Meta materials 


I. INTRODUCTION 


Understanding and using Double Negative (DN) materials— 
negative «(@) and (æ) for some frequencies, in the design of 
devices using electromagnetic waves remains a challenge. 
This topic has become of great interest since these man-made 
materials became popular this decade. Many of us had 
difficulties understanding propagation because we had not 
experienced DN materials before and we had to except the idea 
of a negative phase velocity as a consequence of our 
equations. The concept of a positive phase velocity appears to 
have always been sacrosanct. It’s interesting to note, however, 
that scientists could accept the idea of a phase velocity going 
faster then the speed of light because invariably a real physical 
situation associated with this occurrence could be found. 


Here we derive the propagation equations in DN materials in 
a way that differs from previous derivations. We propose that 
for any real material the energy of the wave must decrease as it 
moves from its source of energy. Energy losses per unit length 
of travel are finite and can be extremely small. Propagation in 
a lossless DN media is found as the mathematical limiting 
solution of an extremely small energy loss per unit length for 
real &(c) and 4(@). We end up with the same conclusions as 
the other authors using nothing more than an algebraic solution 
to Maxwell's equations without the need to even mention 
phase velocity; negative phase velocity is a consequence of our 
solution. 


Why was it necessary to develop this derivation when the 
result was already available? We've gone through the trouble 
because we plan to examine complex multi-dimensional 
structures and do not want to confront potential ambiguities 
about which signs to choose in taking square roots, etc. 
Having a guiding principle to go by is a big help. In this report 
we apply the basic principles to reflection and refraction of an 
electromagnetic waves between a material having positive 
&(o) and (œ), and one having negative &( c) and o). 


In the next section we develop the propagations using the 
principles just mentioned. Use of the equations to solve the 
aforementioned reflection/refraction problem is rendered in 
section 3.Concluding remarks are rendered in section 4 and 
related references are listed in section 5. 


IL. ELECTROMAGNETIC PROAGATION MODEL 


FOR DN MATERIALS 
For a DN material, the Maxwell's equations are 
VeB=0 (1) 
VeD=0 (2) 
vxe=-S 6) 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


r (4) 
ôt 


Using the formulas, 
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(5) 
veo -| 2 Mi ES i+ ey a. 
Cy 0 dS Og; x ) Ox | dy J 
(6) 
we write equations 5 and 6 in component form 
OE, OE,| 0B, (7) 
Oy Oz Ot 
B 
OE, OE, 0B, (8) 
Oz Ox Ot 
E 
OE, E, ôB, (9) 
Ox oy Ot 
OH, 
oH, yay 2 (10) 
ayo ar 
OD, 
Be aly pu (11) 
Oz Ox © ôt 
OH, 
os esr aad (12) 
Ox Oy © 8t 
The constitutive equations are 
D(t) = eE(r) + PO), (13) 
Bt) = u GG) € M) =H) uM) 04 


A dynamic relationship is assumed between the 


polarization P(t) ) and E(t) , and between the magnetization 


M(t)) and H(t) ). An example of dynamic equations for the 


polarization and magnetization is the set provided by Smith 
[16]: 





2P P = - 

= £T, M + OP = £o; E(t) (15) 
"M M B = 

2 £D, x tM = Op E(t) (16) 


The system of equations is solved with the formula (/Ot) = 
joand using the transform pair 


Foe | asco oa (17) 
2 s 


Fa) =f exp jon F cat (18) 


The function F(t) represents either P(t) ; E(t) M (t), 
H(t) , or I(t) . We have 


P(t) = e p E (19) 


DO) -s E-P-e 0 z(0)E- e, (E =6(@)E (20) 
Q1) 
(22) 


£,(0) =1+ y(o) 
E(@) = ££, (@) 


For illustrative purposes, we use the following form for, 
although this is never necessary for calculations. 











qi = ol — ae jor, = us) "E v 
gu (0) = Z's (90) if"; (0) (24) 
0, (@, — O°) 
i LUUD EPOR TO e 
" 0. OU, 
fs sr T (26) 
M) z,H (27) 
B(o) = uy z, CODE Q0) - p) 28) 
ui C9) = (1+ z, C0) (29) 
Ao) = HoH, (@) (30) 
Onp Onp (Ono —O — joy) 
Zu (@)= O.,-@ + jo, - (o2, oy * oT; 
G1) 
Xu (0) = gy (20) jy @) (32) 


As done for, we use an illustrative expression for ry, (c) : 








' ( ) Orp (Ono i a) (33) 
e ee 
4 2 -0 +o Ti 
P Opl y 
Xy (@) (34) 


2 242 272 
(0,,—0 ) +@T,, 


mo 


In Fourier transform space equations 5 and 6 now become 





dE. E, 5 
<2 --jeuh, 35 
z =| jonh, (35) 
OE, OE T 
Lo .LI--ieuH. 36 
[2 2 HS v 
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OE, OF 2 
—-— |--jouH , 37 
| TEES | JouH , (37) 
F. OH.) . E 
is Faas oe J+ joel@)E, (38) 
Oy Oz 
OH, OH | wu - 
=Z- ru =J + jog&(o)E, (39) 
OH, oH E x 
EE Lf =J, + joe(om)E, (40) 
Ox oy É ° 


Now let us consider propagation in the + z direction where 


there are no current sources. Let J =Oforz>0. Thus for 
z>Othe magnitude of the fields must decrease as z 
increases for a lossy medium because all the sources of 
energy lie in the region zx 0. 


For brevity we choose a TEM wave with E, and E , We 


also assume 0/0x 2 0and 0/0y 20. Equations 35-37 and 
38-40 then become 





+ = —jou(o)H , (41) 
OZ 
oH , z 
-—— = jee(o)E, (42) 
Oz 


E(@) = ££, = Eol * Zr (0)) = Ey) (a (0) — jB (o) (43) 





ag (9) — (09 z'g Qo) (44) 

B,(@) =-x" r (0) (45) 

HCO) = Loh, = Holl Xy(@)) = Mo(@y(@) — Pyu C2) (46) 
a(a)=1+ x’, (@)) (47) 

By(@) -1- x", (@) (48) 

We seek solutions of the form 

E, = Eexp(—jkz) 1 (49) 

H, = H exp- jkz) (50) 


and substitute these expressions into equations 41 and 42. We 
get 0/0z = -jk and 


- jkE = - jou(v)H (51) 
jkH = jog(@)E (52) 

A solution occurs only when 
k’ = o €(@) uo) (53) 


All of the equations presented so far have been derived 
before in various forms in most books in electromagnetic wave 
propagation. Solutions for k have been addressed for the 


following three cases: €(@) and u(œ) are both positive, &(c) is 
positive and /(@) is negative, andé(@) is negative and (o) 
is positive. We are now considering the case where both ¢(@) 
and u(@) are both negative. We show that the solution is a 


straightforward extension of the three previous cases to this 
new regime, and needs no additional mathematical or physical 
assumptions other than using the fact that propagation material 
is not loss-free. 


Before going further, a few comments about using the 
foregoing set of equations is appropriate. First, we are dealing 
with casual functions, whose mathematical properties are well 
known. Being causal we need only deal with the w>0 part of 
the complex plane. From equations: 26, 34, 45, and 48 we 
then see that f; (œ) and py (c) are both positive. In addition, 


we also notice that a@,(@) and @,,(@) are even functions of 


c, namely they are functions of @*. These are important 
properties in the analysis. 


Also, even though our analysis uses specific forms for a, (c) , 
Gy (0), Bg (o) and f, (v) ,the aforementioned symmetries 


of the functions are general. 
function. 


They apply to any causal 


Using equations 43 and 46 in equation 53 we get 
b ko (à -jpe May — jBy) 


- ko (QpAy — Be By) — jke (Gp By — s Bx) (54) 


v? 
k = o ey (0) u (o) = px 


(55) 
where c= (u PNTA is the speed of light in vacuum. The 


right hand side of equation 54 is a complex number; therefore 
k will be a complex number. Let us write k as 








k=k (4+ jg) (56) 
where 
: : L2 1/2 
n+ j= (oso, - jPzBu)- jk, (ag Bg — o Ps) 
(57) 
When equation 56 is used in equations 49 and 50 we get 
E, = EexpC- jkz) = EexpC jka) exp(ky&) (58) 
H , = H exp(- jkz) = H expC- jkonz) exp(ko&) (59) 


If the fields are to diminish as z increases, then é must be 
negative. This is how the + or — sign is chosen in equation 
(57). 


For orientation let us consider the traditional case first. This 
is the one where the @ — terms are positive and are much larger 
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than the œ — terms. Equation 57 is then approximated by the E - BG (68) 
equation i 
n+ j£ = aun, - KB, -au B) H = Ha, (69) 
B. BC K = ky PG, (70) 
=+ U2 1.75 Pa , Pe (60) ! . ee eis E : 
(agg) | i a al J In the foregoing equations, a,, a, , and a,are the unit 


i vectors in the x-,y-, and z- directions. Since 7? is 
ee ii Und) ens aio cac qetoM eng tee ne positive the three vectors in equations 68-70 follow the RH 
approximation rule, that is 








1/2 zx = ~ ~ 
(eed Bu l Pr -1 j| Bu i Pr (61) ExH = EHa, (71) 
Og ae 2| ay Qg LES = 
i KxE=k,n© Ea, (72) 
Equation 60 becomes 
Ax K = Hkn”, (73) 
: J 
nt jë ssa bi (Be The RH rule is natural for cases where the phase velocity 
A UH points in the positive z-direction. This is also physically 
comforting since the root mean square Poynting vector, S s 
zat y^ i( as Ba | Oe Be (62) also point in the positive z-direction. Using equations 44 and 
T OEA 1/2 1/2 
2\ ay E 45 we have 


= 1 - =. 1l. a 
If the “+” is used for the right hand side of equation , then S, =~Re(Ex H )-—a, Re(EH )expQk,é O2) (74) 
will be negative and satisfy the conditions of equations 58 and Z 2 











59. This is therefore the correct decision. We then get Using equation 64 we write 
(4) _ 1/2 
ee) (69) Sims = j Re(ExH^)- 5a, Re(EH *) exp(-ky yz) (75) 
1 qi? qi? i 
(+) E H H E 
T 64 + 1 + 
d A ag ap Py up 2 (76) 
1/2 1/2 
Because @;, Œy, Bg, f, are all positive, 7 is positive y? [2 E LT | (77) 
ay Og 


and é“ is negative. Equations 63 and 64 are what we're used 


to seeing. As we shall show, this result leads to the notion of a The last remaining step to calculate S ms iS to compute 


iti h locity in th directi f flow. ~~ » 
p QM C aaa RCM CE Re(EH ) from either equation 51 and 52. Equation 52 gives 


When equations 63 and 64 are inserted into equations 58 and 
: ; . f z _ 0€(@) & 
59 we immediately notice that the phase part, exp(— jkyrz) , is H = = k (78) 





exp(- jk? z) 2 exp - jk, (a a, ) ^ = ed (65) Using equation 43 for ¢(@) and equation 56 for k(@) we 
get 





























ph 
Since the time dependence is exp(joxX) we see that equation f 25:00 - JB, (@)) É (a, — Bz) E E È 
65 defines a frequency dependent phase velocity, kan (0) + J (a) zhi doy pa t i) 8Z, Z 
EUH 2 
(+) d E 
Vu = = (66) 
ph qno (Aray yi (79) 
è " a 
The attenuation is given by exp(k,éz) which is In the foregoing equation we have g =g" with 
1 
: loy all? B all? B (Apay)? -- jy? 
exp(k yg‘ ‘=e 5 | pou ttt (67) pz 2 - pO exp ju 9) (80) 
H E (a, — jpe) 


The foregoing system of equations that we have just Zo = 44/6, is the impedance of free space. P? and y © 
analyzed is called a right-hand (RH) system. In our construct it 
is comprised of the following three vectors: 
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are the amplitude and phase of g=g respectively, and 





(-) 
ph 


T à — j@ 
exp(—jk yy” z) = exp( jk, (ar, |a, |z)? = exp : (92) 
Zs gu is the wave impedance. Since all the parameters V 


Tote Qo: 
are posteo DNASE En S less r TRE eal pantol Since the time dependence is exp(j@t) we see that equation 


g "^is positive. Thus 92 defines a negative frequency dependent phase velocity, 





-2 
inar E| cosy ? 7 C 
Re(EH*)= Bee (81) v =- (93) 
PZ, (a glo, b 

-2 
: E| cosy“ The attenuation is given b 
Sins 7d, | | pog Chr" (82) BRE 

0 








As we see from equations 71-73 the three vectors: points in 


e koz jas? Ba le | Be koZ o 
exp(kģ™ z) = exp 5 EHE Ui =exp 2 y 
LE E [enl |a. 

the positive z direction. E,H , and S,„„form a RH system 





since $ „ms points in the positive direction. Qe 
As before we have 
Now, let us see what happens at the other extreme, when 
both @, and a, are negative. This happens for the set of 
radian frequencies, (c) , that satisfies the conditions i= cT Fe : a = É (95) 
ag (0) =1+ z's (0) «0 (83) o +j) 
Gy (9) =1+ y'y (9) «0 (84) but because of equation 83 this time we need to use the 
: expression 
We write 
iota) gs) E= (e 0) jfi Go) = -e Qa] + jle (96) 
a, (9) = -æy (œ) (86) From equations 89-91 we have 
T E U2 Jj c 
n+ jé z tljarlenl|- bba) + ilara +|@y|Be))") (87) nO + je = {leste ) py | (97) 


Let's first evaluate equation 87 when the // — terms are much Inserting equations 96 and 97 into equation 95 and using 


smaller than the @—terms. This is what we did before and is previous definitions gives 
the usual case of interest. We now get 


. 1/2 1/2 
re dll fie s l ee &J = 


^ les e| 


Again, imposing the condition that ó be negative to ensure 


~- E 





ZO =g Za (99) 


energy decay in the positive z, we now require that we take the 


negative sign of the right hand sign of equation 88. In lieu of e 


Ty 
(e; leu p^ +d 














equations 63 and 64 we now have g?- PUEDE =P exp(jy™) | (100) 
EE 1/2 E = 
n° = -aslen e | nae 
wj2 
1/2 1/2 s E cosy ^ 
go 1| [eel Bu lea Be |__ lo (90) e cle ae are (101) 
2 le; | CA 2 2Z,P 
"t i73 Again we see from equation 101 that the Poynting vector 
y? (E | 2 Bg [en] E | (91) points in the positive z-direction, and therefore, E,H , and 
len | z| CN form a RH system. The fact that this situation is 


accompanied with a negative phase does not violate any law. 


: (o e NE 
Notice that d" andé are both negative: compare Both cases, the double positive and double negative a, and 


equations 20:31 with 70-77 a@,, are described within the same mathematical framework for 


When equations 89-91 are inserted into equations 58 and 59 a plane wave traveling in any direction. 
we immediately notice that the phase part exp(—jk,7z) is 
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Let Q, be the direction of the electric field and Q,,, be the 
direction of the magnetic field. We have 


Q,.Q, -0 (102) 
Q-0,xQ, (103) 
Q, eO, =Q,°Q, -0 (104) 


E = EQ, exp- jky]Q, e r)exp(- jk), er) (105) 
H =F Êy exp(— jkynQ, °F )exp(—jkyjQ, er) (106) 


For double positive materials use 





g? = (aty), (107) 
1/2 172 

ro ol oes | , (108) 
z les 

ge (109) 


and for double negative materials, use 





7? = -aslay (110) 
1/2 1/2 

Zo EE 
H E 

gg (112) 


HI. REFLECTION AND REFRACTION INVOLVING 
DN MATERIAL 


In this section we demonstrate the unique features of 
reflection and refraction involving a DN material. In the three 
cases considered in figures 1 through 3 the electric field is in 
the plane of incidence. The behavior of figure 1 is well 
known; we include this for the reader's orientation. We 
demonstrate the results shown in figures 2 and 3 from the 
solution of the equations. 


In all cases considered in figures 1 through 3 equations 35 to 
40 apply. They are simplified according to the following 
equations beginning with equation 113. . 


Positive € and 
positive u material 





Figure 1 Reflection in air off a positive £ and positive x 
material. 


Negative € and 
negative y material 





Figure 2 Reflection in air off a negative ¢ and negative 
Lu material.. 


Negative € and 
negative u material 





Figure 3 Reflection in air off a negative £ and negative u 
material off an air interference. 
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E OE 2 E ET 5 
OE, E - - jou(o) H , =0 (113) E, 2 —-H, =ZcosgH, (130) 
Oy | Oz : OE 
a ki ~ = 
OE E (114) E, -—-—-H,--ZsingH, (131) 
~——~ |=—jou(o)H Qe l ; 
& à jou(o)H , 
= x u 
OE. OF E Z-zj|-— (132) 
| : -|2-jou(v)H,-0 (15) d 
ox oy Solution for reflection in air off a positive £ and positive 
di of ali u material 
| y à +}. x seco, — (116) 
7 í : From Figure 1, we have the following 
E 5|. joe 20 (117) Q, er = xsin 9, + zcos Q, (133) 
-A T 
e 2 Q, er — xsinQ, + zcos Pp (134) 
OH, 0H,| 0H, jos(a)E, (118) Q, er = xsin Q, + Z COS Qr (135) 
Ox Oy Ox É . 
Using 
eH, : OE, ~ 1 eH . 
= ap = jwe(a) & (119) E, == y (136) 
z k joe Ox 
eH, OE È 
ys z 120 a OH , 
a jose) a (120) É, econ y 137) 
joe Oz 
Inserting equations 119 and 120 into equations 113 to 118 x z 
: niu i: H,,-Hexp ? KQ,er)- H exp @ j(xsin 9, 4 zcos Q) (138) 
EIVes T Vo Vo 
OH, OH - E., =Z,cosg,H exp- Č j(xsi 
_ up J(xsin o,  z cos Q,) (139) 
ze a o" u(o)e(o)H , (121) y 0 Pr p v. 9, Pr 
The solution for equation 121 is E., = Z, cos o, H exp- ^. j(csin 9, + zcos 9,) (140) 





Al, = Aexp|- j(k,x + k.z) |= Aexp- 2(Qe7) (122) a 
) » : ] 
H pr =H exp- j(Q, er) =H exp-— j(xsin 9, +ZCOS Pp) (141) 
Yo Vo 











3 k. ~ 
E, =—H, (123) 
OE ~ oO. . 
" E, r = Zo cos 94H exp- — j(xsin Mp + ZCOS Pp) (142) 
" z v 
E, 2-—-H, (124) ; 
OE ` ~ OR gs 
E 4 E, 4, = Zo cos Øg H exp- — j(xsin Mp + zcos Pz) (143) 
r — xà, + zà, (125) : Vo 
- ~ O.a . a 
Q=sin ga, + cos ga. 126 A, = Hexp- — j(Q; er) = H exp- — j(xsing, + zcosg;) 
ne NOS e a " (144) 
y- zm (127) ~ 


E.r = Zo cos Q4, H exp- j(xsin Q,-zcosQ,) (145) 
0 


a oO. : 
In the foregoing equations “sand ^: are unit vectors in the E.,-Z,cosQ,Hexp-—— j(xsing, +zcosg,) (146) 
=> Vo 
x-and y- directions respectively, X? is the direction of 








; M ec 1 
propagation, 5/97 d , and $54 — Qed. Fora wave Vo = Jestin (147) 
traveling in the positive x- and positive z-direction we have d 
1 
=? sj = 148) 
k, =—sin 128 Vn ( 
geld (128) Jes 
e 
SR e Z, = |^ (149) 
Eo 


Forum for Electromagnetic Research Methods and Application Technologies (FERMAT) 


mem (150) 


In the foregoing expression, H is the magnitude of the 
incident magnetic field, A is the magnitude of the reflection 
coefficient and B is the magnitude of the transmitted wave. 
These quantities are determined by matching boundary 
conditions at z 2 0 


Ë (0,220) H, 0,270) H (0,220) 20 (151) 


H an|- & jGcsin e| + dex) - Z jGcsin e - Bes - Z jasin e| 


0 0 m 


(152) 
E ,(x,2=0)+ E, (220) = E, (2-20) =0 (153) 


Zycoso, H e| - m J(xsin g, ) — Z9 COS QRA a|- x J(xsin r ) = 
Vo Vo 


Z „ cos OB e|- © j(xsin 3 


m 


(154) 


The solution of equations 151 to 154 is available in standard 
electromagnetic texts. We have 


Pr — 9, (155) 
Veok sing, = Je, ui, sin pr (156) 
H+A=B (157) 
Z,cosp,(H —A)=Z,, cosg,B=Z,,cosg,(H +A) (158) 
a Zo COSY, — Z,, COS Py (159) 

Zo COSY, +Z, COS Py 

2Z 

B gener y (160) 





Zo COS 9, +Z,, COS Qs. 


Solution for reflection in air off a negative c and negative 
u material 


The solution for this case is readily found using the 
discussion for the double case combined with the basic 
conclusions deduced in section 2. The basic behavior for 
propagation for a DN material is summarized from equations 
84 to 112 applied to the loss-free model considered here. 
Restating them in the notation of this section we have 


Q, °Q, -O, eQ, -0 (161) 
EE exp(j —Q, er) (162) 
TES j (163) 
H za exp(j — Qï, er) 


Equations 162 and 163 have a negative phase velocity. 


Suppose, for example, that we assume the transmitted wave 
looks like the one for the double positive case just 
considered—that is, figure 1 is valid. What would the solution 
for the system of equations look like? The behavior for the 
incident (equations 138 to 140) and reflected terms (equations 
141 to 143) remains unchanged. However, the equation for 
the transmitted terms is now 


= O o 
H r= Berp cono +zcos @, ) 


m 


(164) 


Bs = Z „ COS, Bexp— j(xsin Q, + ZCco8 pr) 
v 


m 


(165) 


Bs =-Z „sin 9, Bexp Č j(xsin o, + ZCOS 9, ) (166) 
v 


m 


At the boundary, z = 0, the spatial variation of equations 164 
to 166 go like 


Hus > exp j(xsin Pr) 
i v 


m 


(167) 


E,, > exp j(xsing,) 
V 


m 


(168) 


By comparing the foregoing behavior with incident and 
reflected behavior at z=0, which go as 


exp[- (v! v) j(xsin o, )], we see that there is no way a 


transmitted wave traveling in the direction shown in figure 1 
can satisfy the boundary conditions. 


Now let us consider a transmitted wave shown in figure 2. 
For this case, we have 


OQ, er = —xsin 9, + zCos 9, ) (169) 


Again, using the concept of negative phase velocity we now 
have 





H, = Bexp 2 t xsin 9, + zcos øz) (170) 
n v 


m 


Eu = Z m COS p Bexp —- j(-xsin o, +zcosg,) (171) 
v 


m 


(172) 


g 


Es -Z,sngo,;B exp j(-x sin 9, + zcos 9, ) 
V 


m 


Repeating the same procedure as in equations 151 to 154 we 
have 


H en- L j(xsin | + Ae L j(xsin e. = Box - x J(xsin e| 
Vo Vo Vn 
(173) 
oO, " oO. 
Z, cosp, H ex — j(xsin 9, ) — Zo COS 94A ex — j(xsin e)| = 
Yo Vo 


Z, cos QB e- © j(xsin gy ] 
Vn 


(174) 
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The foregoing equations are identical to for negative 
refraction. 


Solution for reflection in a negative ¢ and negative yu 
material off an air interface 


Based on what we learned in the previous cases we can 
readily write down the equations for figure 3 


H,, = H exp join g, + ZCOSQ, ) (175) 


m 


Es = Z „ cos o, H exp —- j(xsin Q, zcosQ,) (176) 
v 


m 


E ,--Z,sing,H exp j(xsin g, + zcosg,) (177) 
v 


Zas 


m 


= QO. . 
H R= ERD d — Z COS Pp) 


m 


(178) 


E = —Z, COS pr Aexp —- j(xsin Dp — ZCOS Pp) (179) 


m 


En =-Z,, sin p,Aexp— j(xsin Pp —ZCOS Pz) (180) 
v 


m 








Hs Bexp 2 JC-xsin 9, + zcos øz) (181) 
à v 


m 


E, = Z, cos 9, Bexp- —- j(-xsin Pr * zcos 9, ) (182) 
V 


m 


E, r =Z, sin 9, Bexp- —- j(-xsin Pr + zcos øz ) (183) 
E V 


m 


Matching boundary conditions at z = 0 gives 


H ez J(xsin e| + sony 2 J(xsin e -B "E j(xsin e| 
Vo Vo 


m 


(184) 
Z,, coso, H e? j(xsin g, ) —Z,, COS QA eq J(xsin Er ) = 
Vo Vo 
Z, cos o, B ed J(xsin ø; ) 
Vin 
(185) 
The solution is 
Pr 7 9, (186) 
JE nhm sin 9, = | Eo Hlo sin Pr (187) 
H+A=B (188) 


Z „coso; (H — A) =Z,cosg,B = Zo coso, (H +A) (189) 


_ Zn COS 9; — Zo COS Pr H 


A (190) 





Z „n COSO; + Zo COS Or 


2Z „ C089, (191) 





Zn COSO; + Zo COS Or 


IV. CONCLUSIONS 


We derive a system of propagation equations and model in a 
Double Negative material in a way that differs from previous 
derivations. Our derivation is based entirely of on the idea that 
real materials always have some loss, and because of this, 
wave energy traveling in a certain direction must always be 
accompanied by a loss of energy in that direction. Additional 
mathematics is not required. Energy losses per unit length of 
travel are finite, and can be extremely small. Propagation in a 
lossless DN media is found as the mathematical limiting 
solution of an extremely small energy loss per unit length. Our 
system of equations and methodology is used to derive the 
reflection and refraction equations between a positive &(«) and 
/{@) and a DN material. These results are in agreement with 
other predictions. 
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